Generating Gowdy cosmological models 
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Using the analogy with stationary axisymmetric solutions, we present a method 
to generate new analytic cosmological solutions of Einstein's equation belonging to 
the class of T 3 Gowdy cosmological models. We show that the solutions can be 
generated from their data at the initial singularity and present the formal general 
solution for arbitrary initial data. We exemplify the method by constructing the 
Kantowski-Sachs cosmological model and a generalization of it that corresponds to 
an unpolarized T 3 Gowdy model. 



I. INTRODUCTION 

Gowdy metrics , which are exact solutions of Einstein's vacuum field equations, 
represent cosmological models with various possible topologies (S 3 , S 1 x S 2 , T 3 ) and are 
spatially compact inhomogeneous spacetimes admitting two commuting spacelike Killing 
vector fields. They are of especial importance for the study of formation of singularities in 
general relativity. Many efforts IH U, 0, 0] have been made in order to derive and analyze 
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this type of solutions, especially by using numerical methods. The physical structure of 
these metrics is sufficiently simple to expect that the singularities can be analyzed in detail, 
but their mathematical structure is sufficiently complicated so that the global dynamical 
behavior is still far from being completely understood. One of the most intriguing questions 
concerns the structure of the curvature singularity that is expected to appear at certain 
spacelike boundary of the associated spacetime. Many studies have been devoted to the so- 
called "asymptotically velocity term dominated" (AVTD) behavior which states that near 
the singularity each point in space is characterized by a different spatiallv homogeneous 
cosmology [8j. The idea of AVTD behavior was originally proposed in Q more than 30 
years ago, but is still being debated. Numerical analysis of the Gowdy models have shown 
that they do become AVTD near the singularity, except at a set of isolated points, where 
there are "spikes" in the behavior of the metric functions. The origin of these spikes is 
investigated in Gowdy metrics have been also analyzed as toy models in quantum 
midisuperspace gravity [3]. 

Numerical methods have been extensively used to investigate Gowdy models, but only 
recently it has been argued that solutions generating techniques can be applied in this case 



H0, 
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12j to generate new solutions and that even a "simple change of coordinates" can 



be applied to rein terp ret certain stationary axisymmetric solutions as S 1 x S 2 Gowdy cos- 
mological models [l3j . The reason why these methods can be used also in this case is due to 
the well-known fact that the solution generating techniques are applicable to any spacetime 
which admits two (commuting) Killing vector fields. In this work, we will concentrate on T 3 
Gowdy cosmological models and will see that a complex coordinate transformation, together 
with a complex change of metric functions, allows us to apply in a straightforward manner 
the well-known solution generating techniques that have been intensively used for stationary 
axisymmetric solutions. 

This paper is organized as follows. In Section |H] we derive the "transformation" that 
relates stationary axisymmetric solutions with Gowdy T 3 models. We will show that due to 
this analogy, the AVTD behavior in Gowdy T 3 is mathematically equivalent to the behavior 
of stationary asixymmetric solutions near the axis of symmetry. In Section IIHI we show 
that Sibgatullin's method for constructing solutions can be applied in the case of Gowdy 
T 3 models and in Section IIVI we present several examples of exact solutions generated by 
using this method. Finally, Section is devoted to the conclusions and some remarks about 
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different possibilities of generalizing the results derived in this work. 

II. STATIONARY AXISYMMETRIC SOLUTIONS AND GOWDY T 3 MODELS 

Consider the line element for stationary axisymmetric spacetimes in the Lewis- 
Papapetrou form 

ds 2 = -e 2 ^(dT + uod<j)) 2 + e- 2 ^[e 2 \dp 2 + dz 2 ) + p 2 d<f) 2 ] , (1) 

where ip, u, and 7 are functions of the nonignorable coordinates p and z. The ignorable 
coordinates T and are associated with the two Killing vector fields rjj = d/dT and r\n = 
d/d<f). The field equations take the form 

1 e 4 ^ 

ip pp + -ip p + il> zz + -^(uj 2 p + co 2 z ) =0 , (2) 
p zp 

lu pp - -LO p + u zz + A(u p ipp + u z ip z ) = , (3) 



e 



'111! 



lp = p{^ p -^ 2 z )-—(u 2 p -u 2 z ) , (4) 

7^ = 2pip p tp z - —e 4 ^u p uj z , (5) 

where the lower indices represent partial derivative with respect to the corresponding coor- 
dinate. 

Consider now the following coordinate transformation (p, t) — > (r, a) and the complex 
change of coordinates (0, z) — > (5, x) defined by 

p = e" T , T = cr, 2 = ix, = (6) 

and introduce the functions P, Q and A by means of the relationships 

*=i(P-r), = », 7 =I(p-£-l). (7) 

Introducing Eqs.tjHJ) and (JJJ) into the line element (JIJ, we obtain 

Let us take r > (what seems reasonable in virtue that the radial coordinate p = e~ T > 0) 
and "compactify" the new coordinates as < x> ^ <$ < 27r (a less reasonable condition since 
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in general — oo < z < +00 and T > 0). The line element (JBJ) with the coordinates r, Xi a an d 
5 in the range given above is known as the line element for Gowdy T 3 cosmological models 
Furthermore, one can verify by direct calculation that the action of the transformations 
(JBJ) and on the field equations 0-© yields exactly the field equations for the Gowdy 
cosmological models which after some algebraic manipulations can be written as a set of two 
second order differential equations for P and Q 

Prr - e~ 2T P xx ~ e 2P (Q 2 T - e-^Ql) = , (9) 

Q TT - e- 2r Q xx + 2(P T Q T - e- 2T P x Q x ) = , (10) 
and two first order differential equations for A 

A T = P 2 T + e^P 2 + e 2P (Q 2 T + e^Q\) , (11) 

X x = 2(P x P T + e 2P Q x Q T ) . (12) 

It should be emphasized that this method for "deriving" the Gowdy line element from 
the stationary axisymmetric one involves real as well as complex transformations at the 
level of coordinates and metric functions. It is, therefore, necessary to demand that the 
resulting metric functions P, Q, and A be real. That means that in general it is not possible 
to take an axisymmetric stationary solution and apply the transformations to obtain a 
Gowdy cosmological model. If the resulting functions are not real, they cannot be physical 
reasonable solutions to the real equations (}9~j)- (fT2"j) . These transformations can be used 
only as a guide to get some insight into the form of the new solutions. In any case, the 
corresponding field equations have to be invoked in order to confirm the correctness of the 
solution. 

A very useful form for analyzing the field equations of spacetimes with two (commuting) 
Killing vector fields is the Ernst representation. In fact, this convenient form for the field 
equations was first proposed for axisymmetric stationary solutions , but since then it has 
been applied in many different configurations If], 0, Q| • Here we will present the Ernst 
representation of the main field equations (jHJ) and (jl(Jj) which is especially adapted to the 
coordinates used here. To this end, let us introduce a new coordinate t and a new function 



R = R(t, x) by means of the equations |10L Il3] 



t = e~ T , R t = te 2P Q x , R x = te 2P Q t . (13) 
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Then, the field equation (jHJ) can be expressed as 



t 2 ( P tt + -P - P xx ) + e- 2P (R 2 t -Rl) = 0, (14) 



whereas Eq. (jl(Jj) for the function Q turns out to be equivalent to the the integrability condi- 
tion R tx = R x t- However, an alternative and convenient equation is obtained by introducing 
Eq.(fT3j) directly into Eq.(fTUj). So we obtain 



te p (Ru + ~Rt - R xx ) - 2[(te p ) t R t - {te p ) x R x } = , (15) 

an equation which of course becomes an identity if the integrability condition R tx = R x t is 
satisfied. We can now introduce the complex Ernst potential E and the complex gradient 
operator D as 

( d d \ 

E = te p + iR 1 and D = ( — , %— 1 , (16) 
which allow us to write the main field equations in the Ernst-like representation 

Re(E) (d 2 E + jDtDE^j - (DE) 2 = . (17) 

It is easy to verify that the field equations (fT^j) and ()15|) can be obtained as the real and 
imaginary part of the Ernst equation (|17|) . respectively. For the sake of completeness, we 
rewrite the system of first order, partial, differential equations (fTT]) and (fE?|) in terms of the 
Ernst potential: 

\t = -Uc + C* + + C-C*_) , (18) 



2 

x x = -- (c + c; - c_ci ) , (19) 



where 



c± = rjW<W-I- (20) 



and the asterisk denotes complex conjugation. 

If the Ernst potential E is known, then it is easy to recover the metric functions P, Q 
and A which enter the line element (jHJ) of Gowdy T 3 cosmological models. In fact, from 
Eq.(|16|) one can algebraically construct the functions P and R. Then the function Q can 
be obtained by solving the system of two first order partial differential equations given in 
(fTHJ) . Notice that the integrability condition of this last system is satisfied by virtue of 
Eq.ffTTj). Finally, the system (fTHj) and (fT^j) for the function A can be solved by quadratures 
since its integrability condition coincides with the Ernst equation (fTTj) . Consequently, all 
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the information about any Gowdy T 3 cosmological model is contained in the corresponding 
Ernst potential. 

One of the most important properties of the Ernst representation (j!7j) is that it is very ap- 
propriate to investigate the symmetries of the field equations. In particular, the symmetries 
of the Ernst equation for stationary axisymmetric spacetimes have been used to develop the 
modern solution generating techniques [19j . like the Backlund method, Belinsky-Zakharov 
inverse scattering method, the Hoenselaers-Kinnersley-Xanthopoulos method, and others 
(for an introductory review and detailed references see I n au these methods it is nec- 

essary to start from a given "seed" solution which has to be specified in the whole spacetime 
(except, perhaps, in the regions where the metric possesses true curvature singularities). 
An alternative approach for exploring the symmetries inherent in the Ernst equation was 



explicitly developed by Sibgatullin 



2l| and consists in constructing exact solutions to the 



Ernst equation from initial data specified only on certain hypersurface (submanifold) of the 
spacetime. For instance, in the case of stationary axisymmetric spacetimes, Sibgatullin's 
method allows one to construct exact solutions from their data on the axis of symmetry. In 
the following sections we will show that Sibgatullin's method can be applied in the case of 
Gowdy cosmological models and will present several examples of its application. 



III. CONSTRUCTING SOLUTIONS FROM AVTD DATA 

As we have mentioned above, an important property of Gowdy cosmological models is 
its AVTD behavior near the initial singularity. In the case of T 3 models it can be shown 
that the singularity is approached in the limit r — > oo. The AVTD behavior implies that 
at the singularity all spatial derivatives of the field equations can be neglected and only the 
temporal behavior is relevant. On the other hand, the transformation © indicates that 
the limit r — > oo is equivalent to the limit p — ^ 0; however, this is true only at the level 
of coordinates and a more detailed analysis is necessary to make sure that this analogy 
is also valid at the level of explicit solutions. To this end, let us consider the system of 
partial differential equations for ip and u given in Eqs.(j2J) and Q. If we neglect the spatial 
dependence on z, which accoring to the transformation (JHJ) is equivalent to the spatial 
dependence on x m Gowdy models, then we obtain a system of differential equations which 
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can be solved by quadratures and yields 

^ = iinK,- +6 y-«)], " = n( ,,4 6 6 V -) +id ' (21) 

where a, 6, c and d are arbitrary real functions of z. Clearly, this solution is meaningless 
when considered as a stationary axisymmetric spacetime. However, if we now follow the 
prescription given in Eqs.(JB} and (JZj) for obtaining Gowdy models, we fill find that solution 
(J2IJ) "corresponds" to the Gowdy model 

P = lnKe- + 6 2 e CT )], Q = ^ + &2) + (22) 

where now a, b, c and d are to be considered as arbitrary real functions of the coordinate 
X- It is straightforward to verify that the expressions given in Eq. (|22j) satisfy the Gowdy 
field equations Q and (|10|) in its "truncated" form, i.e., when the spatial derivatives are 
neglected. The solution (l2*2l is known in the literature as the AVTD solution for Gowdy T 3 

n 

models [6J and dictates the behavior of these models near the singularity r — > oo. Thus, we 
have "derived" the AVTD solution starting from its stationary axisymmetric counterpart. 
This is a further indication that the behavior of Gowdy models at the initial singularity is 
mathematically equivalent to the behavior of stationary axisymmetric solutions at the axis. 
For the sake of completeness we also quote here the value of the function A corresponding 
to the AVTD solution (|22|) that can be obtained by integrating Eq. fjllj) : 

A = A -c 2 lnt, (23) 

where Ao is an additive constant. Furthermore, the corresponding AVTD Ernst potential 
can be obtained by introducing Eq.(|22|) into from Eqs. ()16|) and (fTTj) . Then 

E = a[e" (1+c)T + b 2 e~ {1 - c)r } + iR avtd with R™ td = -2abc . (24) 

If we define 

E(r^oo, X ) = e(x) (25) 

as the Ernst potential at the singularity, we see from Eq. (J24|) that for c G (—1, 1) only the 
imaginary part remains, e(x) = iR avtd . This means that the real part of e(x) is arbitrary 
and since R avtd is given in terms of the real part it is also arbitrary. If c ^ (—1, 1), the 
Ernst potential diverges at the singularity for arbitrary values of the functions a and b. 
In the limiting case c = ±1, the Ernst potential at the singularity is regular, but again 
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no conditions appear for the behavior of the functions a and b. Consequently, the AVTD 
behavior does not impose any conditions on the function e(x)- We will now see that it is 
possible to use this function to construct the corresponding Ernst potential E(t,x)- 

Sibgatullin's method [2]] has been developed to construct exact stationary axisymmetric 
solutions starting from their data on the axis of symmetry. It is based upon the fact that 
the Ernst equation possesses symmetry properties associated with an infinite-dimensional 
Lie group which transforms one solution of this equation into another solution of the same 
equation. This implies remarkable analyticity properties that make it possible to reduce the 
Ernst equation to a system of linear integral equations which can be integrated explicitly 
if initial data is known, for instance, on the axis of symmetry. It is clear that the Ernst- 
like representation (j!7j) possesses similar symmetry properties. On the other hand, we have 
shown that the behavior of stationary axisymmetric solutions near the axis is mathematically 
equivalent to the behavior of Gowdy T 3 cosmological models near the singularity. Thus, it 
should be possible to construct Gowdy cosmological models starting from the value of the 
corresponding Ernst potential at the singularity. It turns out that Sibgatullin's method can 
be generalized in a straightforward manner to include the case of Gowdy models. A detailed 
explanation of the procedure necessary to obtain the system of linear integral equations 



associated with the Ernst equation is given in [2l|. Here we will only quote the main steps 
of the construction. Assume that the value of the Ernst potential is known at the initial 
singularity, i.e. e(x) is given. Then, the Ernst potential can be generated by means of the 
integral equation 

E M = (26) 
where the unknown function has to be found from the singular integral equation 

- 1 mie*(v)+e(0] 



1 (s — k)\/1 — s 2 



ds = 0, (27) 



with the normalization condition 

" x MO 



-ivT 



s 



ds = n, (28) 



where £ = x + ts, r] = x + tK, with s, k 6 [—1, 1]. 

Notice that for this method no condition is imposed on the behavior of e(%). This is in 
accordance with the result obtained above about the AVTD behavior of the Ernst potential 
near the singularity. Once e(x) is given in any desired form, one only has to calculate the 
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integral (|26[) to find the Ernst potential. However, to calculate this integral one first has 
to find the function by means of the singular equation (|27j) and the normalization 
condition (|28|) . In practice, for a given e(£) one has to make a reasonable ansatz for 
such that it allows the existence of solutions for the integral singular equation (|27j) . 



IV. EXAMPLES OF GOWDY T 3 MODELS 



The cases where the Ernst potential at the initial singularity behaves as a rational function 
are relatively easy to analyze. In this section we will present two such examples. Let us 
consider the following simple example of an Ernst potential at the singularity 

e W = , (29) 

where xo is a real constant. The first step of the construction is to find the unknown function 



fj, according to Eqs.(|27jl and (|28j). A reasonable ansatz is again a rational function 21 1 



n = A) + - A- , (30) 
£ — u 

where £i is the root of the equation e(£) + e(£) = (in this case £i = xo) an d A, Al> 

are functions of t and x- To handle the integrals which follow from the singular integral 
equation we use the following standard formulae 

r 1 ds , , 

/-.7T^ = '- (31) 



TT 



i (a + is^VT 3 ^ Va 2 + b 2 

ds 71 



(32) 
(33) 



-i (s - k)(s - 7)Vl - s 2 (k - 7)V7 2 - 1 
where a, 6 and 7 are arbitrary constants. 

Introducing Eq. (j30|) into the normalization condition 1)280. we obtain 

A) + — = 1 , (34) 

r_ 

where r_ = \J (x~ Xo) 2 — whereas the integral singular equation 1)2*7)) leads to 

^ _M±±I=) =0 . (35) 
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The last two equations can be used to find the explicit values of A and A\ which then can 
be replaced in the result of the integration of Eq. fffilj) and yield 

At - 2xqA q 2xo - r+ - r_ 



E(t, x) = -A 



(36) 



r + 2%o + r + + r_ 

where r + = \J (x + Xo) 2 — t 2 - It is easy to check that indeed this is a solution to the Ernst 
equation (fTTj) . Since the resulting Ernst potential is real, the solution corresponds to a 
polarized (Q = 0) Gowdy model. The expression for the metric function P can easily 
be obtained from the definition (|T6*j) and Eq. (j36|) . and the remaining function A can be 
calculated (up to an additive constant) by quadratures from Eq.ffTHj) and (fT9")) : 



A = In 



(r + r_ 



(37) 



t(r + + r_ + 2 Xo )\ 

The physical significance of this solution becomes plausible in a different system of coordi- 
nates which we introduce in two steps. Let us first introduce in the (r, x)-sector of the line 
element (JSJ) coordinates x and y by means of the relationships 



e~ 2T = t 2 



X 2 Q (l-x 2 )(l-y 2 ) 



X = Xoxy 



(38) 



or the inverse transformation law 



r _|_ r _ 



r _ r _ 



x 



2xo 



2xo 



(39) 



so that the metric functions become 

1 — x 

In 



A = In 



(x 2 - y 2 ) 2 



(40) 

X^{l-x*){l-y*){l + x)\ 

and is defined 



\ 1!V /(l-x 2 )(l-y 2 )(l+x) 

The second transformation affects now all the sectors of the line element 
by 

T 

x = 1 , y = cos 6 , a = r , 5 = <p . 

Xo 

Then, after some algebraic manipulations, the metric can be written as 

•2xo 



-ds z = - 



, ) 1 rfT 2 + _ l \ dr 2 + T 2 (d 02 + gin 2 



(41) 



(42) 



T J " V T 

an expression that can immediately be recognized as the Kantowski-Sachs cosmological 



model 



2J]. Thus, we have shown that the Kantowski-Sachs metric can be constructed 



from the value of its Ernst potential at the singularity (}2*9"|) . 
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Consider now the more general case 

ex = — - — — — , 43 

Xo + X + iXi 

where Xo and Xi are rea l constants. The unknown function /i(£) can be sought in the form 

H = Ao + -^ + -^-, (44) 



where £i )2 = ±« = ±y Xo — Xi are the roots of the equation e(£) + e*(£) = 0. Substituting 
Eq.(03J) in the integral equation (f2*Tj) we obtain the system 

-A + ^ + ^ = 0, (45) 

Xo + iXi + ol Xo + *Xi ~ " 

Ai(a + ixi) | A 2 (a-ixi) = Q ^ 



r- (Xo + iXi + a) r + (xo + ixi ~ a) 



where r± = y (x ± a) 2 — t 2 . On the other hand, the normalization condition ()28j) yields 

A + ^ + ^ = l, (47) 

an equation which together with Eqs. ()45|) and (j4*U|) form a closed algebraic system that 
determines the coefficients of the function //: 

a(r + + r_) + 2Xi( r + - r_) 



.4 



o 



a(r + + r_) + ?Xi( r + — r -) + 2 «Xo 



Ai = M«-»Xi)(Xo + % + <») 

a(r + + r_) + ixi( r + - r_) + 2«xo 
. _ r + (a + ixi)(xo + iXi ~ a) 

S±2 



a(r + + r_) + ixi( r + — r -) + 2«Xo 
Finally, we calculate the Ernst potential according to Eq. (J2HJ) and obtain 

771/ . \ A . Ai(xo ~ iXi ~ oc) A 2 (xo-iXi + a) 
E{t,x) = ~A + — — ■ — - + 



r- (Xo + iXi + «) r+ (xo + i'Xi - a) 
2«xo - "(r+ + r_) - ixi{r+ - r_) 



(49) 



2«xo + «( r + + r_) + z%i(r + - r_) 
The calculation of the corresponding metric functions can be carried out as described in the 
last section. When integrating the systems of first order differential equations ()13j) for Q 
and (|18|) and (fT9"j) for A, constants of integration appear which we choose such that a simpler 
representation is obtained in terms of the coordinates used. To write down the final form 
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of the metric functions it is convenient to use the coordinates (x, y) as defined in Eq.(j3HJ) 
with xo replaced by a. Then 

P = ln . xg-aV-xgy 2 j (5Q) 



-x 2 ){l -y 2 )[(xo + ax) 2 + xjy 2 } 



n _ 2xiXo(l - y 2 )(xo + ax) 

V — 2 9 9 2 9 ' V" 31 / 

Xo - a 2 ^ 2 - XIV 

A = ln- Xg(x 2 -y 2 ) 2 (52) 

v/(l-x 2 )(l-y 2 )[(xo + «x) 2 + xfl/ 2 ] 2 
This metric corresponds to an unpolarized generalization of the Kantowski-Sachs cosmolog- 
ical model which is obtained in the limiting case xi = 0- 

It is easy to see that near the singularity (r — > oo) the general Ernst potential 
approaches the corresponding AVTD potential as given in Eq. (|43jl . This could be interpreted 
as an indirect proof of the AVTD behavior of the solution obtained here. A more direct proof 
can be given by analyzing the hyperbolic velocity v = J P 2 + e 2P Q 2 (see Ref.Q). In terms 



of the coordinates x and y, the hyperbolic velocity of the solution (jHOj) and (|5T|) is given as 
a rather cumbersome expression. Nevertheless, it is possible to perform an analysis of its 
behavior by considering the different domains of the coordinates according to the definition 
equation (|39|). One can show that in general < v < 1 what according to Ref.[4] implies 
that the solution is AVTD. 

An important property of the method presented here is that it allows to calculate an 
arbitrary Gowdy T 3 cosmological solution with any degree of accuracy. To this end, let us 
expand in Eq.(J2EJ) the value of the Ernst potential e(£) in Fourier series, 

oo 

e(0 = J2 e fe(^ X) cos(kip) , (53) 

k=0 

where eo = e(x) and we have represented the parameter s as s = cos ip . Then the unknown 
function can be expanded in a similar way 

oo 

MO = Z! X) cos(kip) . (54) 

k=0 



The normalization condition ()28|) can easily be calculated and implies that un = 1. Further- 
more, the general solution of the integral equation (|26J) can be written as |2l| 

Y oo 

E(t, X ) =e(x) + 7;J2 e ^k ■ (55) 
A k=i 
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According to Eq.(|27|). the coefficients [i k have to satisfy the following system of algebraic 
equations 

oo 

Hk(e k +i - e* k+l + e\ k -i\ - e* k _ l{ ) = -2e k . (56) 

1=1 

Thus, once the value of the Ernst potential is given at the initial singularity [e(x) = E(x, t — »■ 
oo))] the general solution of the Ernst equation reduces to an infinite series with coefficients 
satisfying a set of pure algebraic equations. 



V. CONCLUSIONS 

We have shown that it is possible to generate Gowdy T 3 cosmological models starting 
from their data near the initial singularity. To this end, we first show that the Gowdy T 3 
line element can be obtained from the line element of stationary axisymmetric solutions by 
means of complex transformation that involves the metric functions and the coordinates. 
The behavior of stationary axisymmetric solutions at the axis of symmetry is shown to 
be mathematically equivalent to the behavior of Gowdy T 3 models near the singularity. In 
particular, we have derived the AVTD solution from its stationary axisymmetric counterpart. 
We then use the Ernst representation of the field equations and apply Sibgatullin's method 
to the Ernst potential which can be given at the singularity as any arbitrary function of the 
angle coordinate x- I n particular, we have shown that the Kantowski-Sachs cosmological 
model can be derived in this manner by starting from a specific form of the Ernst potential 
in terms of a rational function. We then have found an unpolarized generalization of the 
Kantowski-Sachs cosmological model. This generalization has been obtained in the same way 
as the Kerr metric is obtained from its value at the axis of symmetry by using Sibgatullin's 
method [21] . It is possible to consider more general examples of Ernst potentials at the axis 
in terms of rational functions. It turns out that the system of integral equations (|27j) and 
()28|) forms a closed algebraic system from which the value of the function /i(£) can be found 
and the expression for the Ernst potential can be calculated j^]. This method could also 
be applied in the case of Gowdy cosmological models considered here. 

By expanding the value of the Ernst potential at the singularity in terms of a Fourier 
series, it is possible to write explicitly the general solution for this type of models (including 
the unpolarized case) by using only a recurrence algebraic formula. This is a result that 
could find some application in numerical investigations since it allows us to "control" the 
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accuracy of the analysis by truncating the series at any desired level. 
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